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Abstract 
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1 Introduction 

The main problem of QCD is to understand the structure of hadrons and their properties in 
terms of quarks and gluons. Nucleon charges defined as matrix elements of vector, axial and 
tensor currents between nucleon states contain complete information about quark structure 
of the nucleon. These charges are connected with the leading twist unpolarized q(x), the 
helicity Aq(x) and transversity Sq(x) parton distribution functions (PDFs). The first two 
PDFs have been extensively investigated theoretically and experimentally in many works 
(for instance see [1,2] and references therein as well as [3-5]). There is a big experimental 
problem to measure the transversity of the nucleon because of its chiral odd nature. Only, 
recently the tensor charge 5q(x) was extracted [6] using the data from BELLE [7], HER- 
MES [8] and COMPASS [9] Collaborations. This extraction is based on analysis of the 
measured azimuthal asymmetries in semi-inclusive scattering and those in e + e _ — > h\hiX 
processes. Since Sq(x) is a spin dependent PDF, it is interesting to investigate whether 
there is a "transversity crisis" similar to the case of "spin crisis" in Aq(x). Therefore, 
reliable determination of nucleon tensor charge receives special attention. 

Theoretically, tensor charges of hadrons are studied in different frameworks such as, 
non-relativistic MIT bag model [10], SU(6) quark model [11], quark model with axial 
vector dominance [12], lattice QCD [13], external field [14] and three point versions of QCD 
sum rules [10]. 

In the present work, using the most general form of the nucleon interpolating field, we 
study the tensor form factors of nucleons within light cone QCD sum rules (LCQSR). The 
LCQSR is based on the operator product expansion (OPE) over twist of the operators 
near the light cone, while in the traditional QCD sum rules, the OPE is performed over 
dimensions of the operators. This approach has been widely applied to hadron physics 
(see for example [15]). Note that, the tensor form factors of nucleons up to Q 2 < 1 GeV 2 
(where Q 2 = —q 2 is the Euclidean momentum transfer square) within the SU(3) chiral 
soliton model are studied in [16] (see also [17]). The anomalous tensor form factors are 
studied within the same framework in [18]. These form factors are further studied in lattice 
QCD (see for instance [19]). 

The plan of this paper is as follows. In section 2, we derive sum rules for the tensor 
form factors of the nucleon within LQCSR method. In section 3, we numerically analyze 
the sum rules for the tensor form factors. A comparison of our results on form factors with 
those existing in the literature is also presented in this section. 



2 Light cone sum rules for the nucleon tensor form 
factors 

This section is devoted to derivation of LCQSR for the nucleon tensor form factors. The 
matrix element of the tensor current between initial and final nucleon states is parametrized 
in terms of four form factors as follows [1, 19, 20]: 
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where g M = (p — "P M = + an d <? 2 — — Q 2 - From T-invariance it follows that 
E 1T (Q 2 ) = 0. 

In order to calculate the remaining three tensor form factors within LCQSR, we consider 
the correlation function, 

n^(p,g) = i J d A xe^ x (0 \T{ J N (0) J^(x)}\ N(p)) . (2) 

This correlation function describes transition of the initial nucleon to the final nucleon with 
the help of the tensor current. The most general form of the nucleon interpolating field is 
given as, 
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J N (x) = 2e abc [q Ta {x)CA\q' b {x)j A l 2 q c {x) , (3) 
i=i 

where C is the charge conjugation operator, A\ — I ', A\ = A\ — 75, A\ — t with t being 
an arbitrary parameter and t — — 1 corresponds to the Ioffe current and a, b, c are the color 
indices. The quark flavors are q = u, q = d for the proton and q = d, q = u for the 
neutron. The tensor current is chosen as, 

J^ u = ua^u ± da^d , (4) 

where the upper and lower signs correspond to the isosinglet and isovector cases, respec- 
tively. 

In order to obtain sum rules for the form factors, it is necessary to calculate the cor- 
relation function in terms of quarks and gluons on one side (QCD side), and in terms of 
hadrons on the other side (phenomenological side). These two representations of the cor- 
relation function are then equated. The final step in this method is to apply the Borel 
transformation, which is needed to suppress the higher states and the continuum contribu- 
tions. 

Following this strategy, we start to calculate the phenomenological part. Saturating 
the correlation function with a full set of hadrons carrying the same quantum numbers as 
nucleon and isolating the contributions of the ground state, we get 

n ,,_ (o\J N (o)\N(p>))(N( p i)\UN( P )} 
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where dots stands for contributions of higher states and continuum. The matrix element 
(0 1 7^(0)1 N(p')} entering Eq. (5) is defined as 

(0\j N (0)\N(p')) = \ N u(p) , (6) 

where \n is the residue of the nucleon. Using Eqs. (1), (5) and (6), and performing 
summation over spins of the nucleon, we get, 
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From Eq. (7) we see that there are many structures, and all of them play equal role for 
determination of the tensor form factors of the nucleon. In practical applications, it is more 
useful to work with E T (Q 2 ) = E T {Q 2 ) + 2H T {Q 2 ) rather than E T {Q 2 )- For this reason, we 
choose the structures a^ u , p^q u and p^i for obtaining the sum rules for the form factors 
H T , E T and H T , respectively. 

The correlation function U^ u (p, q) is also calculated in terms of quarks and gluons in 
deep Eucledian domain p 2 = (p — q) 2 << 0. After simple calculations, we get the following 
expression for the correlation function for proton case: 
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Obviously, the correlation function for the neutron case can easily be obtained by making 
the replacement u <->■ d. 

From Eq. (8) it is clear that in order to calculate the correlation function from QCD 
side, we need to know the matrix element, 

Ae abc (Q \u a a (a 1 x)u b T (a 2 x)d a T (a 3 x)\N(p)) , 

where a±, a 2 and a 3 determine the fraction of the nucleon momentum carried by the corre- 
sponding quarks. This matrix element is the main nonperturbative ingredient of the sum 
rules and it is defined in terms of the nucleon distribution amplitudes (DAs). The nucleon 
DAs are studied in detail in [21-23]. 

The light cone expanded expression for the light quark propagator S q (x) is given as, 
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where the mass of the light quarks are neglected, ml = (0.8 ± 0.2) GeV 2 [24] and 
is the gluon field strength tensor. The terms containing G^ v give contributions to four- 
and five-particle distribution functions. These contributions are negligibly small (for more 
detail see [21-23]), and therefore in further analysis, we will neglect these terms. Moreover, 
Borel transformation kills the terms proportional to the quark condensate, and as a result 
only the first term is relevant for our discussion. 

Using the explicit expressions of DAs for the proton and light quark propagators, per- 
forming Fourier transformation and then applying Borel transformation with respect to the 
variable p 2 = (p — q) 2 , which suppresses the contributions of continuum and higher states, 
and choosing the coefficients of the structures a^ v , p^q u and Pp,q v <$, we get the following 
sum rules for the tensor form factors of nucleon: 
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For the form factor E T (Q 2 ) we obtain the following sum rule: 



(11) 



E T {Q 2 ) 



1 



+ 



m N X N \M^J Xo t 2 1 ; ^ ^ 

1 

Q 2 + xgm 2 ^ 
1 Z" 1 dt 

/ T 2 

'xo t 2 



5 



where 



+ 



1 



Q 2 + xlm 2 N 
1 dt 



■e-^ M \l-t)Fl T (x ) 



+ 
+ 
± 
± 



M 2 J X0 t 3 
1 

Q 2 + xlm 2 N 

j_ r 1 ^ 

M 2 L t| ( 



(l-*)^l(*s) 



-/ M2 (l-t)F| T (x ) 



-s(t 2 )/M 2 n _ ,N F 5 



to 
1 



Q 2 + xgm^ 



(l-^I>o) 



1 



Q 2 + rcgm 2 



+ /'f e-^)/- 2 (l-04 T (* 2 ) 



.X() 



/<2 /-A /■!— p _ 

dA J dp J dtiTtfafrl-h-p) , 

(t u p, 1-h-p) , 



Fl T (h) = ~±m z N J dp J dt lP 
"t 2 r l -p 



r 2 + % 



F| T (t 3 ) = 4m 2 N J dp J dt lP \A 2 - V 2 

rl-t2 _ 

F| T (t 2 ) = 8m 2 N / dhTfituh, l-h-t 2 ), 

"l-*3 



(tl,l-tl -P,p) , 



+ Vl M (t 1? l -t! -t 3 ,t 3 ) , 



rl-t-2 

F 7 ~ T (t 2 ) = 8 I dtiTifat^l-ti-h)) , 

•l-*3 



o 



(tl,l-*l-*3,*3) • 



Finally for the form factor H T (Q 2 ) we get the following sum rule: 
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In these expressions, we also use 



F(xi) = x 2 , l-x 1 -x 2 ) , 

^(Xi) = J r (x 1 , 1-Xi- x 3 , x 3 ) , 



s(x, Q 2 ) = (1 — x)m 



(1-x) 



N 



Q 2 , 



X 

where x (s , Q 2 ) is the solution to the equation s(x , Q 2 ) = s . 

The residue Ajv is determined from two-point sum rule. For the general form of the 
interpolating current, it is calculated in [25], whose expression is given as 
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The Borel transformations are implemented by the following subtraction rules [21-23], 
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3 Numerical analysis of the sum rules for the tensor 
form factors of nucleon 

In this section, numerical results of the tensor form factors of nucleon are presented. It 
follows from sum rules for the form factors that the main input parameters are the DAs of 
nucleon, whose explicit expressions and the values of the parameters f^, Xi, A 2 , ff, A™, 
and Vf in the DAs are all given in [21-23]. 

In the numerical analysis, we use two different sets of parameters: 

a) All eight nonperturbative parameters /at, Ai, A 2 , ff, /|, A\ and Vf are estimated 
from QCD sum rules (set 1). 

b) Requiring that all higher conformal spin contributions vanish, fixes five A™, Vf, 
ff, and A%, and the values of the parameters /at, Ai, A 2 are taken from QCD sum rules. 
This set is called asymptotic set or set2. 

The values of all eight nonperturbative parameters (see for example [26]) are presented 
in Table 1, 





Set 1 


Asymptotic set (set2) 


fN 


(5.0 ±0.5) x 10~ 3 


GeV 2 


(5.0 ± 0.5) x 10- 3 GeV 2 


Ai 


(-2.7 ±0.9) x 10 _: 


1 GeV 2 


(-2.7 ±0.9) x 10~ 2 GeV 2 


A 2 


(5.4 ±1.9) x 10" 2 


GeV 2 


(5.4 ± 1.9) x 10~ 2 GeV 2 


Af 


0.38 ±0.15 







V? 


0.23 ±0.03 




1/3 


ft 


0.40 ± 0.05 




1/3 


ft 


0.22 ±0.05 




4/15 


n 


0.07 ±0.05 




1/10 



Table 1: The values of eight input parameters entering the DAs of nucleon. 

The next input parameter of the LCQSR for the tensor form factors is the continuum 
threshold sq. This parameter is determined from the two-point sum rules whose value is 
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in the domain so = (2.25 — 2.50) GeV 2 . The sum rules contain also two extra auxiliary 
parameters, namely Borel parameter M 2 and the parameter t entering the expression of the 
interpolating current for nucleon. Obviously, any physical quantity should be independent 
of these artificial parameters. Therefore, we try to find such regions of M 2 and t, where 
the tensor form factors are insensitive to the variation of these parameters. 

Firstly, we try to obtain the working region of M 2 , where the tensor form factors are 
independent of it, at fixed values of Sq and t. As an example, in Figs. (1) and (2) we present 
the dependence of the tensor form factor H T {Q 2 ) induced by the isoscalar current on M 2 at 
different fixed values of Q 2 and t, and at s = 2.25 GeV 2 and s = 2.50 GeV 2 for sets 1 and 
2, respectively. From these figures, we see that H T (Q 2 ) is practically independent of M 2 at 
fixed values of the parameters Q 2 , sq and t for both sets 1 and 2. Our calculations also show 
that the results are approximately the same for two sets, therefore in further discussion, we 
present the results only for set 1. We perform similar analysis also at Sq = 2.40 GeV 2 and 
observe that the results change maximally about 5%. The upper limit of M 2 is determined 
by requiring that the series of light cone expansion with increasing twist converges, i.e., 
higher twist contributions should be small. Our analysis indeed confirms that the twist-4 
contributions to the sum rules constitute maximally about 8% of the total result when 
M 2 < 2.5 GeV 2 . The lower bound of M 2 is determined by requiring that the contribution 
of the highest power of M 2 is less than, say, 30% of the higher powers of M 2 . Our numerical 
analysis shows that this condition is satisfied when M 2 > 1.0 GeV 2 . Hence, the working 
region of M 2 is decided to be in the interval 1.0 GeV 2 < M 2 < 2.5 GeV 2 . The working 
region of the parameter t is determined in such a way that the tensor form factors are also 
independent of it. Our numerical analysis shows that the form factors are insensitive to 
cos 9 (with t = tan#) when it varies in the region —0.5 < cos 9 < 0.3. 

In Figs. (3)-(5) we present the dependence of the form factors H T (Q 2 ), E T (Q 2 ) and 
H T (Q 2 ) on Q 2 at s = 2.25 GeV 2 , M 2 = 1.2 GeV 2 and fixed values of t, respectively, 
using the central values of all input parameters in set 1 for the isoscalar current. For a 
comparison, we also present the predictions of self consistent chiral soliton model [16] and 
lattice QCD calculations [19, 20] in these figures (note that, chiral soliton model result exists 
only for H T (Q 2 )). 

We see from Fig. (3) that, our results on H T (Q 2 ) are close to the lattice QCD results 
for Q 2 > 2.0 GeV 2 , while the results of two models differ from each other in the region 
1.0 GeV 2 < Q 2 < 2.0 GeV 2 . Our and lattice QCD results differ considerably from the 
predictions of the chiral soliton model. It also follows from these figures that the form 
factors get positive (negative) at negative (positive) values of the parameter t. 

In Figs. (6), (7) and (8), we present the dependence of the form factors H T (Q 2 ) , E T (Q 2 ) 
and H T (Q 2 ) for the isovector current, i.e., for the ua^u — da^ u d current. Our observations 
for set 1 can be summarized as follows: 

• The Q 2 dependence of H T {Q 2 ) is similar to the isoscalar current case, but the values 
are slightly larger compared to the previous case. 

• Similar to the isoscalar case, the form factors Ht{Q 2 ) and Et{Q 2 ) get positive (neg- 
ative) at negative (positive) values of the parameter t. 
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• In contrary to the isoscalar current case, the values of Ht(Q 2 ) are positive (negative) 
for negative (positive) values of t. 

• Our final remark is that the LCQSR results on the form factors can be improved by 
taking into account the a s corrections. 

In conclusion, using the most general form of the nucleon interpolating current, we cal- 
culate the tensor form factors of nucleon within the LCQSR. Our results on these form 
factors are compared with the lattice QCD and chiral soliton model predictions. 

Note added: After completing this work, we become aware of a very recent paper arXiv: 1107.4584 
[hep-ph] [27] in which part of this work is studied. 
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Figure 1: The dependence of the form factor H T of nucleon on M 2 at Q 2 = 1 GeV 2 and 
s = 2.25 GeV 2 , at six different values of t: t = —5; —3; — 1; 1; 3; 5, using the first set of 
DAs for the isoscalar current. 
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Figure 2: The same as in Fig. (1), but at s = 2.5 GeV 2 and using the second set of DAs. 
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Figure 3: The dependence of H T on Q 2 at M 2 = 1.2 Gel/ 2 and s = 2.25 GeV 2 and four 
fixed values of t: t = —5; —3; 3; 5, for the isoscalar current. 
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Figure 4: The same as in Fig. (3), but for the form factor Et{Q 2 )- 
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Figure 5: The same as in Fig. (3), but for the form factor Ht{Q 2 )- 
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Figure 6: The same as in Fig. (3), but for the isovector current. 
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Figure 7: The same as in Fig. (4), but for the isovector current. 
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Figure 8: The same as in Fig. (5), but for the isovector current. 
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